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We derive an effective Hamiltonian for the nonlinear process of parametric down conversion in the 
presence of absorption. Based upon the Green function method for quantizing the electromagnetic 
field, we first set up Heisenberg's equations of motion for a single atom driven by an external 
electric field and in the presence of an absorbing dielectric material. The equations of motion are 
then solved to second order in perturbation theory which, in rotating-wave approximation, yields 
the standard effective interaction Hamiltonian known from free-space nonlinear optics. In a second 
step, we derive the local-field corrected Hamiltonian for an atom embedded in a dielectric host 
medium, i.e. a nonlinear crystal. Here we show that the resulting effective Hamiltonian is found to 
be trilinear in the electric and noise polarization fields, and is thus capable of describing nonlinear 
noise processes. Furthermore, it reduces to the phenomenological nonlinear Hamiltonian for the 
cases where absorption, and hence the noise polarization field, vanishes. 
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I. INTRODUCTION 

Ever since the discovery of second harmonic generation 
by Franken in 1961 [1 , nonlinear optical process have 
been the subject of great interest. Uses of such processes 
cover the full spectrum of possible applications from op- 
tical communications [2] at one end to fundamental tests 
of quantum mechanics [3J Uj at the other. The strongly 
correlated photons that are created in these processes are 
regularly used in many quantum cryptographic protocols 
[5ji6j, and in the areas of quantum information processing 
and quantum computing [7j. As a result, these process 
have been the subject of much study (for a necessarily 
incomplete selection, see e.g. Refs. (SHU])- 

The fundamental theory that describes the interaction 
of light and matter is quantum electrodynamics (QED). 
This theory has proven to be highly successful over the 
past sixty years, and has accurately described a disparate 
range of physical phenomena over a wide range of ener- 
gies, from scattering of charged particles in high energy 
colliders to the low energy dynamics of atoms in elec- 
tromagnetic fields. The theory which, in its microscopic 
form, describes the interaction of electromagnetic fields 
with charged particles, predicts the appearance of non- 
linear processes when high intensity fields interact with 
certain types of matter. However, owing to the complex- 
ities associated with the microscopic structure of matter, 
the calculation of such properties are highly involved and 
often neglect some of the more complicated features. For 
example, in such calculations absorption is almost always 
neglected. As a result the standard approach to nonlinear 
optical processes is, in the main, phenomenological and, 
although it provides a good approximation at high in- 
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tensities [12], does not necessarily hold for all situations. 
Furthermore, there are several circumstances where non- 
linear absorption is thought to be a critical factor [T3] 
and hence the standard approach is insufficient in these 
cases. 



It is certainly the case that absorption plays an im- 
portant role in many physical processes. Recently, a 
method for consistently quantizing the electromagnetic 
field in absorbing linear electric and magnetic materials 
has been developed (for reviews, see e.g. Refs. [El US]). 
There have been some attempts to extend this theory to 
nonlinear materials [T71[T5] , however, as yet a full theory 
has proven to be elusive. 

In this article we present an extension of the linear 
quantum theory of light in absorbing media to nonlinear 
processes. In the following we will consider, as an ex- 
ample, the second order process of parametric down con- 
version, where an input (pump) photon is converted by a 
nonlinear medium to give two output photons (signal and 
idler) whose frequencies sum to that of the input photon. 
We will begin, in Sec. [ll] by briefly reviewing the quan- 
tization scheme for linear absorbing dielectric materials. 
In Sec. |III[ we consider the interaction of photons with a 
single atom to second order in perturbation theory and 
derive an effective interaction Hamiltonian for this pro- 
cess. In Sec. |IV[ in order to find the nonlinear response 
for a bulk material, we apply local field corrections to 
the Green functions of the interacting electric fields by 
considering the effect of placing the interacting atom in 
a cavity within a bulk material. This method produces 
an effective interaction Hamiltonian for the second or- 
der nonlinear process that includes the sought nonlinear 
noise processes. Concluding remarks are given in Sec. \V\ 
Some useful expressions and lengthy derivations can be 
found in the Appendices. 
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II. ELECTROMAGNETIC FIELD 
QUANTIZATION IN LINEAR MEDIA 

Before we outline the theory of electromagnetic field 
quantization in nonlinear dielectric media we shall briefly 
review the quantization of the electromagnetic field in a 
linearly responding medium [15[ I16j . We begin with the 
classical Maxwell equations in frequency space. In the 
absence of free currents and charges these equations take 
the form 



with 



V-B(r,c^)=0, 
V X E(r, uj) - za;B(r, w) = 0, 
V-D(r,c.) = 0, 
V X B(r, w) + za;^oD(r, w) = 0, 



D(r,w) =£oE(r,w) + P(r,cj) 



(1) 

(2) 
(3) 
(4) 



(5) 



Although magnetoelectric media can be treated along the 
same lines, for simplicity we assume that there are no 
magnetic responses present. For spatially local isotropic 
media, the general form of the linear polarization field is 

P(r, uj) = 60X^1) (r, c.)E(r, c.) + PN(r, c.). (6) 

The first term in Eq. Q is the linear response of the 
medium to an external electric field with a linear suscep- 
tibility x*^^-'(r,w). The second term is the linear noise 
polarization field which describes Langevin noise that is 
associated with absorption, and is required for the theory 
to be consistent with the fluctuation-dissipation theorem. 
As a result the frequency components of the electric field 
obey the inhomogeneous Helmholtz equation 

^2 

VxVxE(r,cj) ^£(r,cj)E(r,w) = [i^V^iv.uj), (7) 

where £(r, cj) = \ + ^^^{v,lS) is the complex permittivity 
of the medium. This equation can be formally solved 
using the Green tensor for the Helmholtz operator 

E(r, uS) ^ cjVo y" d^s G(r, s, cj) • Pn(s, uj). (8) 

The Green tensor G(r, s, iS) solves the Helmholtz equa- 
tion with a point source 

^2 

V X V X G(r,s,w) ^e(r,a;)G(r, s,w) = <5(r-s). (9) 

Quantization is then performed by relating the noise 
polarization field to a set of bosonic field operators 



PN(r,c.)=z^^£"(r,c.)f(r,c.) (10) 
and imposing canonical commutation relations for them, 
f(r,a;),f^(s,w')l = ^(r - s)5(w - w'). (11) 



Thus the frequency components of the quantized electric 
field can be written as 

E(r,w)=iW^^ / d?s^e"{s,u:)G{T,s,u) ■iis,^) 

(12) 

and the total field operator reads 



E(r) / dwE(r,cj) + h.c. 



(13) 



The bosonic operators f(s,a;) and f^^(s,a;) describe col- 
lective excitations of the electromagnetic field and the 
absorbing dielectric material and can be viewed as the 
generalization of the free space photonic mode operators 
to arbitrary media. The bilinear Hamiltonian 



oo 

Hf= J d^r J duhuji'f{r,uj) -{{i 



w), 



(14) 



generates the time-dependent Maxwell equations from 
Heisenberg's equations of motion for the electromagnetic 
field operators. 

This quantization scheme has already been successfully 
applied to a wide range of linear media including mag- 
netic and magnetodielectric materials and has been used 
to study many physical effects including spontaneous re- 
laxation rates, atom surface interactions and various cav- 
ity QED processes [16]. 



III. EQUATIONS OF MOTION OF THE 
LIGHT-ATOM SYSTEM 



Starting with this theory, it is possible to study nonlin- 
ear interactions of light with an atom in the presence of 
dielectric bodies. Here, one solves the coupled equations 
of motion recursively to obtain an expansion in powers 
of the electric field operator. Each of these higher-order 
terms corresponds to a specific nonlinear process. From 
these interaction terms effective nonlinear Hamiltonians 
can be derived which characterise each of these processes. 
In this section we will look at the derivation of the effec- 
tive Hamiltonian for the second order nonlinear process 
of parametric down conversion, where an input (pump) 
photon with frequency ujp is converted by a nonlinear 
crystal to give two output photons (a signal photon with 
frequency uj^ and an idler photon with frequency w^) such 
that ujp ^ ujs + uJi. 

The interaction between light and a single atom can 
be described using the multipolar coupling in the dipole 
approximation. Here the electromagnetic field couples 
linearly to the dipole moment of the atom. The multipo- 
lar coupling Hamiltonian can be written as 



H — Hy + H ^ + H\, 



(15) 
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with 

i 

Hint = d ■ E(r^) 



'3- ' diu' 



'J 

X v/£"(s,c^')GpA(rA, s, A(s, c^') + h.c. (16) 



and Hp from Eq. (14) above. The above Hamiltonian 



has been written in component form with the Greek in- 
dices running over the three Cartesian coordinates. To 
these indices the summation convention apphes. The 
atomic Hamiltonian is the sum of the projectors onto 
the (undisturbed) energy levels of the free atom located 
at r^, a-ii = is the projector onto the ith eigenstate 
with energy fkUi and aij = is the atomic flip oper- 

ator between the iih. and jth atomic energy state. Note 
that we do not apply the index summation convention 
over the (Latin) atomic state indices. 

From the Hamiltonian ( 15 ) we obtain Heisenberg's 



equations of motion for the atomic and bosonic field op- 
erators as 



^ I d\s I duo 



X |[5Ajfc(rA,s,cj)(7ifc - gx,k2{'rA,s,uj)&kj] /a(s,w) 

(18) 

where LOij = uji~ ujj are the atomic transition frequencies 
and the coupling constants g\,ij{r,s,Ltj) are defined by 



gx,ijir,s,uj) 



(19) 



The differential equations (17) and (18) completely de- 



scribe the dynamics of the coupled light-atom system. 
As we are primarily interested in the dynamics of the ra- 
diation field, the next step will be to remove the atomic 
degrees of freedom. This is done by formally solving for 
the atomic operators and resubstituting the result into 
the equation for the field operators. The result will be 
a single dynamic equation for the radiation field in the 
presence of the atom. The formal solution of Eq. (18) is 



dt' 



dh 



duj e*""' 



,(*-*') 



X |[gAjfc(rA,s,w)(7ife(t') ~ gx,ki{rA,s,uj)akjit')] fx{s,uj) 

+ bAjfe(rA,s,w)<Tjfc(t') -.gI,,fc(rA,s,a;)<Tfcj(i')] /l(s,t^)} 

(20) 

where (Jij(O) is the atomic operator at t — 0. This is 
a recursive expression for aij. The full solution can be 



obtained by substituting Eq. (20) back into itself. The 



result is an infinite series in increasing powers of the field 
operators fx{s,LL!) and /|(s,aj). Each of these higher 
order terms corresponds to a specific nonlinear process 
that can occur when the radiation field interacts with the 
atom. We are interested in parametric down conversion 
and hence are interested in the term that is quadratic in 
the field operators. In fact, since the parametric down 
conversion process creates two photons we will be inter- 
ested in terms quadratic in field creation operators. It 
transpires that terms containing the annihilation opera- 
tor will average to zero when we apply the rotating wave 
approximation. In the light of this, for notational clar- 
ity, we will henceforth not display in detail terms that 
contain either /a(s,w) and/or /^(s',a;'). Resubstituting 
Eq. (20) twice and then inserting it back into Eq. (17) 
gives 



ft (r, io") = iLo"fl (r, - * E Su., {VA, r, l 

ijkp 



t t' 




^glfe,(r^,s,c.) (a.fe(O)e--*' -*e--(*'-*") [^^^^.(r^, s', c.')^.p(t") - S^^.plr^, s', c.')^pfe(i")] /^(s',^')) 
gl,,,(r^,s,..) (a,, (0)e-^^-*' -ze--^^- (*'-*") [5;,,,(r^,s',^')<^fcp(t") " 5;;op(rA,s',a.')ap,(t")] ^(s',^.')) }/l(s,^) 



terms containing /a(s, w) and/or /^(s', oj') 



(21) 
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It is worth momentarily digressing from the derivation 
to consider the consistency of Eq. (21 ) with the coupled 
Eqs. (17) and (18). Here we have formally solved the 



atomic equations of motion and substituted the result 
into the equation of motion for the electromagnetic field, 
thereby describing the effect of the atom on the field. In 
order to complete the analysis, one also has to study the 
backreaction of the field on the atom. As a result of this 
backreaction, the bare atomic transition frequencies ujij 
become complex-valued quantities, gaining a line width 



Fjj and a level shift Sujij 



Uij + Sujij 



(22) 



The derivation of this result and the explicit expressions 
for the the line shifts SuJij and the line widths Tij, in 
terms of the dyadic Green function, is well documented 
in the literature (see, e.g. Refs. [THl HO] for reviews). 
In reality this feature will not be critical in the follow- 
ing derivation, but it is important when considering the 
causality properties of the result (see Appendix[B]). How- 
ever, for the consistency of Eq. (21 1 with the original 



equations of motion, the modified transition frequencies 



ujij in Eq. ( 22 ) must be used in place of their bare coun- 
terparts. 



We now return to Eq. (21). In order to solve the in- 



tegrals a number of approximations need to be made. 
Firstly we write the bosonic field operators as the prod- 
uct of a rapidly oscillating function and a slowly varying 
envelope function. 



(23) 



Secondly we assume that the radiation is off-resonant 
with any of the atomic transitions and hence the fre- 



quency of the radiation field and those associated with 
the atomic transitions are significantly different. 

We now apply the rotating-wave approximation. The 
Ihs of Eq. (21) evolves at a frequency w". Terms on 



the rhs of Eq. (21 1 contain contributions from the two 



field modes evolving at uj and uj' and contributions from 
the atomic operators. We are interested in the process 
of parametric down conversion where the frequencies of 
the incoming and outgoing photons combine such that 
uj" = oj' + u!. Terms on the rhs whose frequencies of 
evolution deviate significantly from the above condition 
will oscillate rapidly in comparison to the resonant terms. 
Hence, over long time periods these non-resonant terms 
average to zero. Thus, we keep all terms on the rhs of 
Eq. ( 21 ) that satisfy uj" — oj' + uj. As a result we neglect 



terms which contain the annihilation operator fx{s,uj). 

Furthermore, in order to obey w" = w'-l-w, none of the 
(far off-resonant) atomic transition frequencies uJij can 
appear on the rhs. Hence, we drop all off-diagonal atomic 
operators and retain only those terms that contain diag- 
onal atomic projection operators an. Lastly terms with 
aij (0) represent the free (undriven) motion of the atom 
in the background field and hence are not of interest here. 

The physical motivation for these approximations 
comes from the nature of the parametric down conver- 
sion process itself. The condition w" = uj' + uj is a state- 
ment of energy conservation. Hence there is no energy 
available to drive atomic transitions and thus the atom 
must stay in its initial state. Therefore, the off diago- 
nal atomic operators, which describe atomic transitions, 
cannot contribute. 

After applying these simplifications and permutating 
some indices we find 



t t' 

fl{Y,Uj")^iJ'fl{v,J')+l J dt' J dt" J J dPs' j duj J duj'Y^^r. 

"^3^ 



^)5;,H(rA, s', uj')g,A^A, r, uj")] /^(s, ^)/t (s', uj'). 



r 



(24) 



We would now like to perform the time integrals. Note 
here that the coupling constants g\^ij{v,s,uj) are not 
functions of time and hence can be taken out of the in- 
tegral. Since the slowly varying envelope of the field op- 
erator /j^(s,a;) is approximately constant over the time 
periods of interest, it can also be taken out of the integral. 



Lastly the atomic operators an are the projection opera- 
tors on to the energy eigenstates of the atomic Hamilto- 
nian and hence stationary under evolution by the atomic 
Hamiltonian. Thus, the first term in Eq. (24) integrates 
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to 



t t' 







Using Eq. ( 23 ) we can recombine the rapidly varying part 



of the bosonic operators with the slowly varying envelope 
to recover the full time dependent operator 



(w' - - UJkj) (w' - CJife)(w + UJ' - LOij) 



x/l(s,c.)/t(s',c.') 



1 



I , w ^ , r/l(s,c.)/t(s',c.'), (26) 

where we have neglected the rapidly oscillating terms 
since these will again average to zero over long time peri- 
ods. Integrating the other terms in Eq. (24) in a similar 
way gives 



/t(r,c.") = *c.'7t(r,w") d's d\s' / du; / rfc.' 



X i^A^,(r^;s,s',r;c^,c^',c^")/l(s,c^)4(s',c^'), (27) 

where we defined the nonlinear coupling tensor operator 
K),f,^{rA;s,s',r;uj,uj' ,uj") as 

ijk 



(w' -UJik){uJ+Uj' -UJtj) 

9*x,ijii^A, s, a;)ff* fc,(rA, s', uj')g^^kj{TrA, r, uj") 

{uj' - OJki)iuJ + W' - Wfcj) 

91m(^a, s, uj)g*^,^{rA, s', uj')g^,kjirA, r, uj") 

{uj' - UJij){uj +UJ' - UJkj) 

al.jkjT^A, s, uj)g*^j^.^{rA, s', uj')g^^j,{rA, r, uj") 

{uj' - UJki){uJ + UJ' - UJji) 



(28) 



It is now straightforward to write down an effective 
interaction Hamiltonian that, via Heisenberg's equations 
of motion, generates the correct dynamical equation for 
the bosonic field operators: 



= -h I d^r I d^s I d^s' I duj / duj' / duj" 



X 7^AM.(rA; s, s', r; J ,uj") fl{s, uj)fl{s', uj')l{r, uj") 

+ h.c.. (29) 

It is evident that the dynamical evolution of the atomic 
quantities is frozen out in this approximation. The non- 
linear coupling tensor operator depends solely on the pro- 
jection operators an on to the atomic eigenstates. Since 



their evolution is now static, (Ti 



them by their expectation values p[^\ Hence, the effec- 



0, we can replace 



five interaction Hamiltonian ( 29 ) becomes a functional of 



the dynamical variables of the quantized electromagnetic 
field alone. 

Although the interaction part of the Hamiltonian in 



(25) 

p9|) correctly describes this process at a microscopic level 



in terms of bosonic operators, it is the macroscopic de- 
scription, in terms of electric fields, that is of practical 
interest. The response of an atom to an applied electric 
field is described in terms of a single atom susceptibil- 
ity or polarizability. Although an essential part of the 
macroscopic description, the polarizability is, in fact, a 
function of the microscopic properties of the atom. Hence 
we can use this to relate the microscopic effective Hamil- 
tonian to an equivalent macroscopic effective Hamilto- 
nian. It can be shown (see Appendices |A] and |b| that 
the causal second order nonlinear polarizability can be 
written as 



xS-,(w,a;') 



1 



(0) 



ijk 



d(x ,ik dj3 ,j i ^7 , kj 



da .jidj^^ik dj _ i^j .kjdp^ik dj j" i 



{uj' - UJki){uJ +UJ' - UJkj) {uj' - UJki){uj + UJ' - UJji) 

(30) 

One should also note that in using this form of the po- 
larizability requires the condition uj" = uj + uj' . Using 
this we can write Eq. (28 1 in terms of the second order 
nonlinear polarizability 



//N _ « (2) r . , ./^ / ^^0 



Kxf,^{rA;s,s',r;uj,uj',uj") = jXapti^,^') 

-y/e"{s,uj)e"{s',uj')e"{v,uj") 



X G*^y^{YA,S,Uj)G*p^{vA,s' ,Uj')G^y{vA,V,Uj"). 

(31) 

Hence the interaction term of the effective Hamiltonian 
becomes 



HI 



cfF 



d^r / d^s / d^s' / duj / duj' / duj" 



xxS,(^.^') 

X Ga\{'rA-,&,uj)G*p^{TA,s' ,uj')G^i,{ya,y,uj") 

X fl{s,uj)fl{s',uj')f,{v,uj")+h.c.. (32) 

Finally, we can combine the Green functions with the 
various factors to re-form electric field operators 



,6 £.2 



-y/e"{s,uj)e"{s' ,uj')e"{r,uj") 



duj I duj' 



duj"x^^l^{co,co') 



X Ei{rA,to)El{rA,Lu')E^{rA,uj") + h.c. . (33) 
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The Hamiltonian ( 33 ) describes a second order nonlinear 



interaction between the quantized electromagnetic field 
and a single atom in free space, possibly near (but out- 
side) a dielectric body. In order to describe the situation 
in which the nonlinearly responding atom is located in- 
side a dielectric or even part of the dielectric medium 
itself, an additional ingredient is necessary. 



IV. LOCAL FIELD CORRECTIONS 

Previously, we have considered the nonlinear interac- 
tion of an electric field and a single atom in free space. 
In this case the applied fields act directly on the atom, 
and thus the local field at the atom &°'^{rA) is equal to 
the applied field Ea{^A)- In the case where the interact- 
ing atom is part of a larger body, the electric field at the 
position the atom is different from the applied external 
field. The surrounding material modifies the applied field 
such that 



(34) 



The local field correction method involves calculating the 
prefactor £[£(cj)] so that the local interaction can be re- 
lated to the applied fields. This is a common technique in 
linear optics and has even been applied to nonlinear pro- 
cesses [21] . There are a number of ways to perform these 
corrections. Here we shall consider the real cavity model, 
which was first discussed in the framework of quantum 
optics in Ref. j22| and in its present form in Ref. [53]. 
This technique is well known and has already been used 
to calculate a number of atomic properties such as mod- 
ified spontaneous decay rates (24j and one and two atom 
van der Waals interactions [25] . 

In this model the interacting atom is placed inside an 
empty spherical cavity of radius Rc, which itself is em- 
bedded in the host medium. The local field correction is 
performed by replacing the Green function found in the 
expansion of the electric field with that of the spherical 
cavity. The Green function for the spherical cavity can be 
found by considering wave propagation from the cavity 
centred at ta to a point r located in the host medium. 
This is similar to the Onsager model [IS] for local 
field corrections (Appendix [C|), a technique that is more 
common in classical nonlinear optics, where the correc- 
tions to the classical fields are calculated using similar 
concepts. The Green function method, however, is more 
general as it takes into account the absorptive properties 
of the surrounding material whereas the Onsager model 
does not. 

Consider an atom at the centre of an empty spherical 
cavity of radius Rc embedded in an infinite medium of 
permittivity e{uj) such that 

^('^'^)-\e(..) if |r-r^|>i?, ^^^> 

with Rc on the order of the interatomic distance. The 
Green function for the spherical cavity can be split up 




Rfc,r,co) 



FIG. 1: The interacting atom is placed at the centre of a 
cavity of radius Rc, embedded in an infinite homogeneous 
medium. The scattering part of the Green function can be 
split into a contribution Rap{rA, r, i^) related to reflection off 
the cavity wall and a contribution Tai3{rA,r,(jj) related to 
transmission through the cavity wall. 



into two parts 

Gai3{rA,r, uj) = C^pivA, r, u) + Gf^(rA, r, w), (36) 

where G^^(r^,r,a;) is the part which describes trans- 
mission within the cavity medium (i.e. free space) and 
Gf^(rA,r,a;) is the part which describes scattering off 
the cavity wall. Furthermore, we use the decomposition 

G^p{rA,r,uj) = Rai}{rA,r,uj) + Tafi{rA,r,uj) (37) 

where i?Q,^(r^, r, w) is the contribution from reflection 
off the cavity wall and TQ^(r^,r,a;) is the contribution 
from transmission through the cavity wall. We shall 
consider a coarse-grained model with the characteristic 
length scale much greater than Rc, the interatomic dis- 
tance. Thus the individual atoms cannot be resolved and 
so the medium can be considered to be a uniform and 
continuous. Furthermore, since the cavity is of radius 
Rc there is only one resolvable point within the cavity, 
the location of the atom, ta- As a result there can be 
no propagation within the cavity and thus G^^(r^,r,w) 
can be neglected. 

The reflective part of the Green function is [211 |25l [28] 



Rc,^{rA,r,uj) = '^C{uj4TTRlS{rA-r)5c^f, (38) 

DTT 6 



with the Mie reflection coefficient 



G(c.) = 



eicj)h['\z) zoj\''{zo) -Jl^'izo) zhY'iz) 



(39) 

where Zq = luRc/c and z = e{ui)ijjRc/ c and j''i\z) and 
h^i\z) are, respectively, the first spherical Bessel and 
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Hankel functions of the first kind, 

Z \Z Z'' J 

(40) 

As Rc is small compared to the optical wavelengths as- 
sociated with the process (i.e. Rc <^ c/uj) we can expand 
C{uj) in powers of ujRc/c 



C{uj) = 3 



e{io) - 1 

;2e(w) + 1] iuj^Rl 
9 f [ejuj) - l][4e{Lj) + 1] \ c 
5 \ [2e{Lo) + 1]2 J iwi?, 

. r, £(a^)^^^(^) _ 1 _^ c"^'^^^ 



[2e(w) + 1]2 



Substituting these expressions into i?Q^(r^,r,w) and 
taking the cavity radius to zero {ujRc/c — ?> 0) gives 

R.,{rA, r,.) = l 0^S^{-A - r)S^,. (42) 
The transmission part of the Green function is [211 [25l 



Tc,p{rA,r,uj) = D{io)G^p{rA,r,uj) (43) 
with the Mie transmission coefficient 



Diu:) 



j['\zo) \zoh['\zo)] -h['\zo) IzoA'Hzo) 



j«(zo) \zh^^\z)\ -e{i,)h['\z) \zoj['Hzo) 



(44) 

Here, G^^(r^,r,u;) is the bulk Green function for the 
infinitely extended medium without cavities. Expanding 
D{u!) in powers of uiRc/c gives 



3e(w) 



O 



ujR, 



(45) 



2e(w) + 1 

which, in the limit coRc/c — > 0, leads to 

T„^(r^,r,..) ^ -^£M_G^^(r^,r,c.) (46) 

for the transmission part of the Green function. Collect- 
ing all the results yields 

^ ^ 2 eiuj) - 1 

G„MrA,r,..) = -^-^-^_5(r^ -r)5„, 

for the local-field corrected Green function including the 
cavity. One should note that the J- function in Eq. (47) 



is the contribution from the reflective part of the scat- 
tering Green function. The intrinsic singularity present 
in the bulk Green function is still contained within 



Here we have used the Green function for a cavity em- 
bedded in an infinite homogeneous media to locally field 
correct the electric field at the location of the atom. How- 
ever, it was shown in Ref. that this result can be eas- 
ily generalized to cavities embedded in arbitrary media. 
To obtain the expression for arbitrary media one merely 
needs to replace the bulk Green function and the homoge- 
neous permittivity with the appropriate Green function 
and inhomogeneous permittivity for the new geometry. 
Hence the following results are, with the right substitu- 
tions, correct for any system. 

We can now apply these local-field corrections to the 
nonlinear Hamiltonian. Inserting the local-field corrected 
Green function ( 47 ) into the effective nonlinear Hamilto- 



(41) nian ( |33[ ) gives 



duj I doj' / doj' 



TT J 



(2) / 



,6 £.2 



Ve"(s,w)£"(s',cj')e"(r,cj") 



C*{u)^6{rA - r)<5„A + i?*(c<.)G^J(r^, r, w) 



C{oj")4j^6{va - v")&,, + ^(c.")G^,(r^, r", J') 



x/l(r,c.)/t(r',c.')/.(r",c.") + h.c. 



where we abbreviated 
2 e(w) - 1 



3£(a;) 



^ ' 3 2e{uj) + 1 ' ^ ' 2e{uj) 

X2) 



(48) 



(49) 



Recall that ^^^^(w,^;') is the polarizability for an iso- 



lated atom in free space and that cj" = uj 



Inside 



a material body, the interpretation of the nonlinear po- 
larizability has changed. We can thus define a local-field 
corrected polarizability (i.e. the polarizability of an atom 
embedded in an extended medium) as 



X^X^u^.u^') - D*HD*{JW')x'^l^{u:.^'). (50) 

Performing the integrations over the (5-functions and re- 
combining the various factors to form electric and linear 
noise polarization fields gives 



-nint - ^0 



duj j duj' j duj"x"^l^{uj,uj') 

Ei{vA.uj)+C[e*{u)]Pl^{vA,u) 

El{vA,uj') + C[e*{u')]Plf,{vAM) 
E^iVA, LO") + C[e{Lu")]P^,^{vA, LO") 



h.c. 



where 



2 e(a;)- 
9eo e{uj) 



(51) 



(52) 
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is the local correction factor for the noise polarization 
field. 

The interaction Hamiltonian ( |51[ ) is trilinear in the 
electric field and the (local-field corrected) noise polar- 
ization field with the strength of the interaction charac- 
terised by the (local-field corrected) second order polar- 

izability Xai^^{uj,uj'). This can be viewed as the combi- 
nation of several different types of interaction processes. 
The term trilinear in the electric field corresponds to the 
sought parametric down-conversion process. All other 
terms describe various nonlinear interactions between the 
electric field and noise polarization fields. They corre- 
spond to the absorption of one or both of the outgoing 
photons by the medium, the production of one or two 
outgoing photons from an excited noise field (e.g. by 
thermal excitation) and a pure nonlinear noise field in- 
teraction. These extra interaction terms are features of 
the absorptive properties of the medium and are in effect 
corrections to lower-order (nonlinear) processes. 

On taking the limit of vanishing absorption the noise 
polarization fields vanish identically [recall Eq. (10)]. 



Hence, the interaction Hamiltonian (51) reduces to the 
standard interaction Hamiltonian associated with para- 
metric down conversion in non-absorbing media. It 
should be noted that computing higher-order nonlinear 
processes will similarly lead to corrections to second order 
processes. However, in view of the applied rotating-wave 
approximation, their effect remains negligible. 



noise polarization fields. 

The local-field corrections have two effects: to re- 
place the free-space second order nonlinear polarizability 

(2) (2) 

Xai3-yi^^^') by its bulk modified counterpart xl^^^{oj,u}'), 
and to introduce additional contributions to the effective 
interaction Hamiltonian. These additional terms are re- 
lated to nonlinear absorption processes inside the bulk 
material, similar to those found in the purely macroscopic 
approach pursued in |181 119| . In contrast to the macro- 
scopic picture, in our present microscopic derivation we 
have gained a better understanding of the origins of the 
additional contributions to the effective Hamiltonian. 

In the limit of vanishing absorption, where the noise 
polarization field disappears, one recovers the standard 
second order effective interaction Hamiltonian as used in 
classical nonlinear optics. In the generic situation when 
absorption cannot be disregarded, the effective Hamilto- 
nian (51 ) will be the starting point for subsequent inves- 



tigations into the role of absorption on the generation of 
down-converted photons and their propagation through 
nonlinear media. 
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V. SUMMARY 

Beginning with the total Hamiltonian for the medium- 
assisted electromagnetic field interaction with a single 
A^- level atom (15), we derived Heisenberg's equations of 



motion for the dynamical variables of both the atom and 
the medium-assisted field. By integrating out the atomic 
degrees of freedom and then focusing on one particular 
term in the nonlinear expansion of the field variables, we 
have found an effective equation of motion for one par- 
ticular optical process. One can think of this equation 
of motion as the dynamical equation that is generated 
from an effective Hamiltonian which describes only the 
process of interest. This effective Hamiltonian does not 
refer to the A^-level atom, consisting only of the dynam- 
ical variables relating to the medium-assisted field. All 
information about the effect of the atom is now contained 
in the coupling constant for the interaction terms. 

Here, as a specific example, by studying the equa- 
tions of motion to second order in perturbation theory 
we have derived an effective Hamiltonian for the second 
order nonlinear process of parametric down conversion 
in the presence of an absorbing host material. By ap- 
plying local-field corrections to the Hamiltonian we have 
moved from considering the interaction of light with a 
single atom to considering the interaction of light with a 
bulk material. We found that the Hamiltonian can be ex- 
pressed in terms of various products of electric and linear 



Appendix A: Second order polarizability 

Here we derive the second order susceptibility of a sin- 
gle atom (also known as the second order polarizabil- 
ity) in terms of the microscopic properties of the atom. 
A semi-classical approach with classical radiation and a 
quantized atom is used. This approach follows closely 
the derivation in Ref. 

In general media the susceptibility x is a tensorial func- 
tion that relates the polarization of the medium to the 
strength of the applied electric field. In the time domain, 
for the case of second order processes, the polarization is 
related to the applied electric field by 

(r,t) = 

oo oo 

eojdrj dT'x^^l^{T,T')E^{v,t-T)Ep{v,t-T'). (Al) 

T 

(2) 

Here Xafi-y ^he second order susceptibility. If the 
medium consists of a single atom, Eq. ( Al I still ap- 



plies with the left hand side giving the polarization of 
the single atom. The susceptibility that appears in the 
single atom version of Eq. (Al) is often referred to as 



the second order polarizability of the atom. 

Consider a single atom in state p which evolves un- 
der the perturbed Hamiltonian H — Ha + Hint- Since 
the operator p can be written as a sum over the projec- 
tion operators onto the energy eigenstates of the atomic 
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Hamiltonian, the state is time stationary under evolution 
of the unperturbed atomic Hamiltonian Hf^. Hence the 
Heisenberg equation of motion for the atom is 



dm 



1 



(A2) 



dt ih ^ 

This can be formally solved to give the recursive relation 
t 

m = dti[p{h),HnAti)]+p^"'> (A3) 

to 

where is a constant of integration and is equal to 
the state of the system at i = to- Equation (A3) can be 
resubstituted into itself to give a series solution for p{t), 

X [[•••[p("\ffi„t(t„)]--- , Hint{t2)h Hintih)] +p 



(0) 



Given that we now have an expression for the quantum 
state of the atom, we can now write down the polarization 
of the atom, which is defined as the expectation value of 
the dipole moment operator 



P^{t) = {d-,) = Tr p{t)d-y{t) 



(A4) 



As p{t) is an infinite expansion so P-y{t) will be an infinite 
expansion. Here second order processes are of interest so 
only the second order term is considered. 



t ti 

dh /"dt2Tr{[[p("),ffi„t(t2)],Hint(ii)] dm] 



io io 



(A5) 



In the dipole approximation, H-mtit) is given by 

H-nm ^ -L{t)E^{t). (A6) 



Hence (A5| becomes 
1 



t ti 

dti J dt2 

to 



X Tr 



{ , -d^ {t2)Ep {t2%-L{ti)Em)\ dj it) } 



{my 



t ti 

dti j dt2 

to to 



X Tr 



(A7) 



where the cyclicity of the trace and the classical nature of 
the radiation have been used to arrive at the result. We 
now make the substitutions ti = t — ti and t2 ^ t — T2 
and take — >■ — cxi (and hence ti,2 — >= +oo). Hence 



Eq. (A7| becomes 



dri / c?r; 



X Tr 



oo oo 

[[d^{t)J^{t - Tl)]Jp{t - T2)] } 

y. E^{t-Ti)Ep{t-T2). (A8) 



The time dependence of quantum operators can be writ- 
ten as 0[t) = W{t)Otj(t) where U is the unitary oper- 
ator U = e-^f^t/h associated with the Hamiltonian H. 
Hence, we can write 



d^{t~T)^U\t)d^{-T)U{t). 



(A9) 



Noting that the dipole moment operator commutes with 
Hint and that p*^"-' commutes with Hx and again using 
the cyclicity of the trace, Eq. ( A8 1 becomes 



oo oo 



Pf(t) = ^/^n/rf.. 





X Tr 



[[dm)Ja{~ri)]Jp{-r2)]} 

y. E^{t-Ti)E0{t~T2). (AlO) 



Comparing Eq. (AlO I with Eq. (All one can see that 



the second order nonlinear polarizability is given by 



(2) ( \ 



1 1 



Tr- 



[K(0),d„(-ri)],d^(-r2)]} 



(All) 



In the frequency domain, the Fourier transforms of the 
electric field and the polarization field are substituted 
into Eq. (AlO), and we find that 



oo oo 



1 



dTi I dT2e 



X Tr \^dm),L{-ri)]Jp{.-r2)\]E^{u)Em'). 

(A12) 



and hence 



1 1 



oo oo 



dri / dT2e ^"'^e 



X Tr- 



{pW [K(0),J„(-Ti)],ciV-r2)]} 



(A13) 
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Note here that, by performing the Fourier transform, we 
have imposed the condition w" = w + w'. The trace of 
the combination of dipole moment operators can be com- 
puted by inserting identity operators in terms of atomic 
energy eigenstates \ei) . The matrix elements of the dipole 
moment operators are given by 

{e,\L{t)\ej) = {e,\U\t)d^tj{t)\ej) 

= {e,\e'"^'/^do,e~'"^'/^\ej) = d„,,je'"-^"*. (AM) 

Here is the energy eigenvalue associated with the eigen- 
state Note also that p^*^^ is diagonal in the eigenbasis 
of the atomic Hamiltonian Hx^ 



ij "u- 



(A15) 



Using the relations (A14), (A15) and the completeness 



relation for the atomic eigenstates, the second order sus- 
ceptibility becomes 



Xa/3T,(W,W') ^ I dn I dT2 



1 1 



ijk 



da J k , ki ^j.ij ^ 



da^kidj^ ^ijd'^ jl^E ^ 6 ^ ^ ^ 



(A16) 

It is important to note that, in order for the expression 
to be consistent with the atom-field equations of motion, 
the transition frequencies uJij must be complex variables 



iVi 



(A17) 



As in Eq. (22), ujij is the bare atomic transition fre- 



quency, 6ujij is the level shift and Tij is the transition 
linewidth. This results in a factor of e"^"^ in the inte- 
grand which leads to convergence at the upper limit of 
the integral. Integrating both time integrals leads to 



(2) / M 1 1 



(0) 



ijk 



doL^ij dp ^ ki dj j' k 



{UJ' ~ U}ik){uJ +UJ' - CJjj) (CJ' - UJik){u} +U}' - UJjk) 



da.,kidp ^ijd^ jli 



doL .jkdp^ij d^ ^ ^ 2 



[uj' - i^ji){i^ + a;' - ujjk) (cj' - ^3i){^ + cj' - uJki) 

(A18) 

which is the expression for the polarizability of a single 
atom in frequency space. 



Appendix B: Causality and the Kramers-Kronig 
Relations 



The linear polarization field is a linear response to the 
applied electric field. As with any response theory the 



magnitude of the reaction is described by the response 
function. In the case of the linear polarization field the 
response function is the linear susceptibility 



oo 



(Bl) 



By causality xi^jC''") must vanish for r < 0; the polar- 
ization field at time t cannot depend on electric fields at 
times greater than t. Thus 



xS(r) = e(r)x2(r), 



(B2) 



where 0(t) is the Heaviside step function. Fourier trans- 
forming both sides of the equation gives 



1 
27ri 



du 



(B3) 



The function \/{uj — uj') has to be seen in its distribu- 
tional sense and, by using Sochotzki's formula, can be 
decomposed into l/(a; — w') = V / [uj ~ lo') + mS / {w — uj') . 
Here V denotes the principal part. This results in 



,xS(^') 



dw 



(B4) 



Decomposing Eq. (B4) into its real and imaginary parts 
gives 



Re 



Im 



xS(^) 



TT 



duj'- 



Im 



xS(-') 



xS(^) 



TT 



duj' - 



ReJxS(-2 



(B5) 



(B6) 



These are the linear Kramers-Kronig relations. A causal 
linear response function must satisfy these relations. For 
more information about the linear Kramers-Kronig rela- 
tions the reader is referred to Ref. (30j . 

One can derive a second order nonlinear version of the 
Kramers-Kronig relations in a similar way. The second 
order nonlinear polarization is given by 

OO OO 

eojdrj dr' x^^l^{T.T')Ep{v,t - T)E^{r,t - t'), (B7) 

T 

(2') 

where XQ^-y('''i '''') is the second order response function. 
Owing to causality, Xap-yi'^i'''') must satisfy 

xLl(T,r') = e(r)0(r')xil(r,r'). (B8) 
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By Fourier transforming Eq. (B8), we find that 

i2) 



37r2 J 

— OO 



StTZ 7 [UJ — Uj) 6711 J {OJ' — UJ') 

— OO — OO 



,(2) 



(B9) 



Splitting up Eq. ( B9 1 into its real and imaginary parts 



gives the nonlinear Kramers-Kronig relations of the form 

(2) / '\ 



Re 



V 

37r2 



dCjduj 



Re 




Co') 




-uj){Cj' - 





v_ 

3^ 



Im 



tl) — w 



v_ 

Stt 



Im 



, (BIO) 



Im 



37r2 



Im 



dwd£b' 



[Cj — uj){uj — LO') 



Stt 



Re 



7^ 



Re 



W' - OJ' 



(Bll) 



As in the case of linear response functions, all causal sec- 
ond order response function must satisfy these relations. 



The polarizability of an atom is a response function 
and hence must satisfy the relevant Kramers-Kronig re- 
lations. The full expression for the second order nonlinear 
polarizability in absorbing media is 



■ijk 



1 1 





— ujik — Suik 


- iTik){uj + Lo' 

da . ij df^ . ki ^7 J k 




- *r,j) 




iOik — Scuik - 


iTik){uj + uj' - 

da^kidp ^ij d'y^jk 


ujjk — Sujjk - 


- i^jk) 




OOji — ^^ji — 


iTji){u! + oj' — 

dajkdp^ijd'yki 


LOjk — Sujjk - 


- i^jk) 


(^'- 


Ulji — Sujji — 


iTji){uj + Lo' — 


ujki — SuJki — 


i^ki). 



(B12) 

which must obey Eq. (B9|. Note that the response func- 



tion is constructed from four terms of the general form 



Tahrf(w,tj') 



A 



(w - UJab - *rafc)(w + UJ' - LOad - i'^'ad) ' 

(B13) 

where ^ is a constant. Applying Eq. (B9) to (B13[) gives 



Tabd[to,Lu ) = -T-^ / dujduj — — — + — / duj—z ^ ^ "-UJ 



371-2 7_oo ((Ii-a;)(w'-a;') 37r {uj - uj) Stt J^oo (w' - w') ' 

■p roc. ^ 

rQbd(w,a;') = -— ^ / dCodu}' 



37^2 (a) - Uj){Ld' - - LOab - i'^ab)i^ + id' - UJad " i^ad) 

iV f°° _ A 



duj- 



y_oo -Uj){u}- LUab - iT ab){i^ + w' - LOad - ad) 

iv r A 



371- '^^ (w' - uj'){u - UJab - «rQb)(a; + Co' - uJad - iT^ad) ' ^^^^^ 



The integrals in (B14) can be solved by residue calcu- 
lus to show that the rhs is indeed identically Ta{,d(cL>,w'). 
Thus terms of this type, and hence the second order non- 
linear polarizability, obey the relevant Kramers-Kronig 
relations. Therefore, the second order nonlinear polariz- 



ability is a causal response function. 
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Appendix C: The Onsager Model 

The Onsager model for local field corrections, as used 
in classical nonlinear optics, is the classical variant of 
the real cavity model [HI [53] and involves considering a 
point charge at the centre of a empty spherical cavity em- 
bedded within a dielectric medium. One then considers 
two cases: firstly, the field in the empty cavity when the 
dielectric is subjected to an external electric field and, 
secondly, the field in the cavity as a result of the po- 
larization the point charge induces in the surrounding 
dielectric. By considering the Maxwell equations for the 
field on the boundary of the cavity, one can derive re- 
lations between the applied and local fields electric and 
polarization fields, 



E 



loc " 

2 

3eo 



2e{uj) + 1 



E, 



2e{uj) + 1 



(CI) 
(C2) 



One should also note that since Ptot = P + Pn , the noise 
polarization field can be assumed to be corrected in the 
same way as the reactive polarization field. 



This method results in the same local field correction 
factors as the real cavity model used for quantum sys- 
tems. For a more detailed description of the method 
readers are referred to Onsager's original paper [3B]. A 
good summary of the model can also be found in Ref. [27] . 
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